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Abstract 

We show that the exact energy eigenvalues and eigenfunctions of the 
Schrodinger equation for charged particles moving in certain class of non- 
central potentials can be easily calculated analytically in a simple and elegant 
manner by using Nikiforov and Uvarov (NU) method. We discuss the gen- 
eralized Coulomb and harmonic oscillator systems. We study the Hartmann 
Coulomb and the ring-shaped and compound Coulomb plus Aharanov-Bohm 
potentials as special cases. The results are in exact agreement with other 
methods. 
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I. INTRODUCTION 



A large number of important physical problems require solving the non-relativistic and rel- 
ativistic wave equations for charged particles moving in central and non-central potentials 
to determine the energy eigenvalues and the corresponding eigenf unctions. It is known that 
for very limited potentials, the Schrodinger equation is exactly solvable. In general, one 
has to resort to numerical techniques or approximation schemes. For many of the quantum 
mechanical systems, most popular approximation methods such as shifted 1/N expansion, 
WKB method, perturbation theory, path integral solution, supersymmetric quantum me- 
chanics and the idea of shape invariance are widely used for this purpose. But some of these 
methods have drawbacks in application. Although some methods give simple relations for 
the eigenvalues, however, they give very complicated relations for the eigenf unctions. 

The study of exact solutions of the Schrodinger equation for a class of non-central po- 
tentials with a vector potential and a non-central scalar potential is of considerable interest 
in quantum chemistry [1-10]. In recent years, numerous studies [11-15] have been made in 
analyzing the bound states of an electron in a Coulomb field with simultaneous presence 
of Aharanov-Bohm (AB) [16] field, and/or a magnetic Dirac monopole [17], and Aharanov- 
Bohm plus oscillator (ABO) systems. In most of these studies, the eigenvalues and eigen- 
functions are obtained via seperation of variables in spherical or other orthogonal curvilinear 
coordinate systems. Recently, the path integral for particles moving in non-central poten- 
tials is evaluated to derive the energy spectrum of this system analytically [18]. Further, the 
idea of supersymmetry and shape invariance is also used to obtain exact solutions of such 
non-central but seperable potentials [19]. The radial and angular pieces of the Laplacian of 
the Schrodinger equation can be treated by using the idea of shape invariance within the 
same framework. 

Very recently, an alternative method called as Nikiforov-Uvarov (NU) method, which 
received much interest, has been introduced for solving Schrodinger equation [20] for some 
well known analytically solvable potentials [21-26]. Similarly, Dirac, Klein-Gordon equations 
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for Coulomb potential [27,28], generalized Hulthen [29,30], usual and modified Woods-Saxon 
potentials are also solved by using the NU method [26,31-33]. 

This paper is structured as follows: In Section II, we briefly introduce the basic concepts 
of the Nikiforov-Uvarov (NU) method. In Section III, we have obtained the eigenvalues and 
eigenfunctions of the generalized Coulomb system and generalized oscillatory potentials in 
the Schrodinger equation using the NU method. In Section IV, we give some concluding 
remarks. 



II. BASIC CONCEPTS OF THE METHOD 

The solution of the Schrodinger-like second-order differential equations plays a vital role 
in studying many important problems of theoretical physics. In this regard, the Nikiforov 
and Uvarov (NU) method can be used to solve these types of equations with an appropriate 
coordinate transformation s = s(r) [20]: 

"^"^"^ + ^/-^'^ + ^"^"^"^ = °' 
where a{s) and cr(s) are the polynomials with at most of second- degree, and f(s) is a first- 
degree polynomial. The special orthogonal polynomials [20] reduce Eq.(l) to a simple form 
by employing ipn{s) — (pn{s)yn{s), and choosing an appropriate function 0„(s). Consequently, 
Eq.(l) can be reduced into an equation of the following hypergeometric type: 

<^(«)2/:(«) + r{s)y',{s) + Ay„(s) = 0, (2) 

where r(s) = f(s) -|- 27r(s) (its derivative must be negative) and A is a constant given in the 
form 

A = A„ = -nr'is) - ~ ^K " (s) , n = 0, 1, 2, ... (3) 

It is worthwhile to note that A or A^ are obtained from a particular solution of the form 
y{s) = yn{s) which is a polynomial of degree n. Further, yn{s) is the hypergeometric-type 
function whose polynomial solutions are given by Rodrigues relation 
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where Bn is the normahzation constant and the weight function p{s) must satisfy the con- 
dition [20] 

d t(s) 

—w{s) = ——w{s), w{s) = a{s)p{s). (5) 
as a[s) 

In order to determine the weight function given in Eq.(5), we must obtain the following 
polynomial: 



In principle, the expression under the square root sign in Eq.(6) can be arranged as the 
square of a polynomial. This is possible only if its discriminant is zero. In this case, an 
equation for k is obtained. After solving this equation, the obtained values of k are included 
in the NU method and here there is a relationship between A and A; by A; = A — 7r'(s). After 
this point an appropriate (f)n{s) can be extracted from the condition 

0'(g) ^ ^ /yN 



III. APPLICATIONS 

Non-central potentials are normally not discussed in most text books on quantum me- 
chanics. This is presumably because most of them are not analytically solvable. However, 
it is worth noting that there is a class of non-central potentials in three dimensions (3D) for 
which the Schrodinger equation is seperable. This section deals with two systems involving a 
generalized form of such potentials: a one known in quantum chemistry as the Hartmann sys- 
tem and another one of much use in quantum chemistry and nuclear physics. Both systems 
correspond to ring-shaped potentials; namely, the Coulomb-kepler and the isotropic har- 
monic oscillator system in 3D. Three-dimensional potentials that are singular along curves 
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have received a great deal of attention in recent years. In particlar, the Coulombic ring- 
shaped potential [1-4] revived in quantum chemistry by Hartmann and coworkers [5,6] and 
the oscillatory ring-shaped potential [7-9] systematically studied by Quesne [10] have been 
investigated from a quantum mechanical point of view by using various approaches. As 
ring-shaped systems, they may play an important role in all situations. For example, the 
Coulombic- Kepler system is of interest for ring-shaped molecules like cyclic polyenes [5,6]. 
Further, the harmonic oscillator system is of potential use in the study of super-deformed 
nuclei. 

In the following we solve the Schrodinger equation in 3D for its eigenvalues and eigen- 
functions with non-central but seperable potentials of much intrest. 



A. The generalized Coulomb system 

The Coulombic ring-shaped (Hartmann) potential [5,6] in spherical coordinates is 

T./ n^ ^ B Ccosd 
r sm 9 sm 9 

where A = —Ze^, B and C are real constants. The potential in (8) introduced by Makarov 
et al. [34] and its importance lies on the fact that compound Coulomb plus Aharanov- 
Bohm potential [16] and Hartmann ring-shaped potential, originally proposed as model for 
the benzene molecule [5] are mathematically linked to this potential. In fact the energy 
spectrum for these two potentials can be obtained directly by considering these as special 
case of the general non-central seperable potential. 

Our aim is to derive the energy spectrum for a moving charged particle in the presence 
of a potential given by Eq.(8) analytically in a very simple way. The Schrodinger equation 
in spherical polar coordinates for a particle in the presence of a potential V{r, 9) is 

«Hi a r.aA 1 / 1 a r i a^ 



2ii 



dr \ dr J r^ \sm9 d9 \ 89 J sm^ 9 d^p^ ^ 

+ V^(r,^)-£;j*„,^(r,^,(^) = 0, (9) 
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where /i — ^^^^^ being the reduced mass and '^nim{fi f) being the total wave function 
defined by 



*nim(r, 9, if) = i?(r)e(^)e±^™^ R{r) = g{r)/r, m = 0, 1, 2, 



(10) 



The Schrodinger equation in Eq.(9) can be reduced to the two ordinary differential equations, 
one for an electron in a Coulomb-like field: 



d'^R(r) 2 dR 

— H h 

dr^ r dr 



2m 



En 



A 
r 



R{r) = 0, 



(11) 



and one allows us to obtain the standard properties of the spherical harmonics given by the 
angular part: 



d^e{e) cose dQ{e) 
de^ ^ sine de 



+ 



[m 



B + Ccos( 



sin'e 



(12) 



We shall first solve the radial part of the Schrodinger equation. Hence, by using the trans- 
formation s — r, Eq.(ll) can be transformed into the form given by Eq.(l) as 



R"{s) + -R'{s) + 
s 



-<S2 _ _ ^ ;l) 



R{s) = 0, 



(13) 



where 



9 2mEn , „ s ^9 2mA 
nl = > (£;„ < 0), (5^ = 



n 



2 ' 



(14) 



representing the bound-state case. Further, defining the following associated polynomials: 



f (s) = 2, a(s) = s, a(s) = - (3^s - 1(1 + 1), 

and using Eq.(6), 7r(s) is found in four possible values as [26,30-32] 

KnS + 1, -P^ + {21 + l)Kn, 

-{k^s + 1 + 1), ki = -P^ + (21 + l)/t„, 

7r(s) = < 

KnS-l - 1, k2 = -P"^ - {21 + 1)k„, 
-KnS + 1, -P^ - {21 + l)Kn, 



(15) 



(16) 
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r s 



(17) 



where the roots are determined by means of the same procedures as in [26,30-32]. We also 
find 

2(«;„s + / + l), A;i = + (2i + 
-2(k„s + 0, A;i = + (2Z + 
2Ks-/), A;2 = - (2/ + 
-2{KnS - / - 1), A;2 - - (2/ + 
Imposing the condition r'is) < 0, one gets the following physical particular solution: for the 
case k2 = — — {21 + we obtain A„ = 2n/t„ and A = — + 2(/ + With A = A^, 
the energy eigenvalues obtained from Eq.(ll) for the Coulomb potential (^4 = — Ze^) are 

-2ii{ZeY 



E, 



n,l 



r, n= 0,1,2, 



(18) 



[2% + / + l)]'^ 

which is the desired bound-state satisfying the boundary conditions [35] . To find eigenfunc- 
tions, we first determine the weight function which satisfies the conditions in Eq.(15) and 
Eq.(17), using Eq.(5) we get 



p{r) = e-''^"V+\ (0 < r < oo). 
Substituting Eq.(19) into Eq.(4), we obtain the first part of the wavefunction 

y„(s) = C„Lf+i)(2«;„r), 



(19) 



(20) 



where L„(r) being the Laguerre Polynomials and C„ being the normalization factor. Using 
Eq.(7), we find the other part of the wavefunction as 



(21) 



Finally, for bound states, the radial part of the wavefunction (finite for all r) has the form: 



Rn{r) = iV„r'e-''»^Lf +i)(2«:,r), (0 < r < oo). 



(22) 



where Nn being a new normalization constant. It is clear from Eq.(22) that the radial wave 
function does vanish at the origin (r = 0) [35]. Hence, Eq.(18) and Eq.(22) are true bound 
state solutions [35]. 



Repeating for the angular part by using the simple transformation s — cos 9, one can 
rewrite Eq.(12) in the following form 

^ - + (T^ [-'(' + - - m = 0. (23) 

Further, defining the following associated polynomials 

f(s) = -2s, (7(s) = l-s^ a(s) ^ -1(1 + l)s'^ -Cs-rri^ - B + 1(1 + 1), (24) 

and using Eq.(6), 7r(s) can be found in four possible values as 

y/A^s + ^/A^, ki ^ l{l + 1) - Ai, 
-^/As + ^/A^, h^l{l + l)-Ai, 

^/A^S + ^/A^, k2 = l{l + 1) - ^2, 

-^/Ms + ^/A[, k2 ^ l{l + I) - A2, 

where A, = ^^^±eW^±^^^^ and A^ = (!!!!±W(^WE^^ being strictly positive 
constants. The roots are determined by means of the same previous procedures as before. 
We also find ris) as 



7r(s) 



(25) 



t{s) = < 



2 

-2 
2 

-2 



\y/A-l)s + ^/A 

/]4l+l)s- 
'(^/T2 -l)s + ^/T, 
+ 1) s - ^/A 

Imposing the condition r'(s) < 0, we also get the following two particular solutions: 

i) For the case ki — l{l + 1) — Ai, one obtains X — l{l + 1) — Ai — \J~Ax and A~ 
n(n + 1) + 2n\/]47, which consequently give the following quantum number: 



A;i = /(/ + 1) - Ai, 

A;2 = Z(Z+1)-A2, 

k2^l{l + \)- A2. 



(26) 



I ^n + JAi, 



(27) 



where n is a non-negative integer. 

Following similar procedures in finding the radial part of the wave function, we also find 
the angular part as 
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(28) 

with \/A^ = —n + /; —n — I — 1. 
ii) For the case k2 = l{l + 1) — ^2, one obtains \ — l{l + 1) — A2 — and A~ = 
ri(ri + 1) + 2n\/~A2i which consequently give the following quantum number: 

l^n+ yi^, (29) 

where n is a non-negative integer. Moreover, the angular part of the wavefunction reads 
e(^) = (1 - cos^)(^-^)/2(l + cos^)(^+^)/'Pi/^-^'^+^V), ^(<^) = e±^^^, 

(30) 

with = -n + l\-n - I - I. For the case (S = C = 0, = 0, ^2 = m^), Eq.(30) 
becomes 

e(^) = sin"* ePiT'"'\cose), $((^) = e^^'"'^, (Sl) 

with m = —n + I; —n — I — 1. 

Therefore, our final energy levels and eigenfunctions for a real bound charged particle in 
a Coulombic field potential plus a combination of non-central potentials given by Eq.(12) 

are 



E. 



n,n,m 



2r 



n 



1/2- 



2 ' 



(32) 



and 



'^nnm{r,0,'p) = N^r^s'^^' L^^'+'\2Knr){l - cos^)(^-^)/^ 



x(l + cos^)(^+^/2pi/^-^'^+^V)e^'^^ 



(33) 



respectively, with the angular quantum number satisfying the condition given by Eq.(29). 
These results are consistent with that in Ref.[19]. Further, for C = 0, it agrees with the 
corresponding energy spectrum already been calculated by the path integral solution of the 
system [18] and using KS transformation [36,37] 



1. Hartmann ring-shaped potential 



The Hartmann ring-shaped potential in question is of the form [38] 

V(r,S) = r,o^{^-^-^\e„, (34) 

\ r sm 6 J 

1 2 4 

where o^o — and Eq — are the Bohr's radius and H-atom ground state energy, 

respectively. The dimensionless parameters rj and a are positive and real which range from 
about 1 to 10 in molecular applications. This potential had been proposed by Hartmann as 
a model for molecules like benzene [5] and can be obtained from the non-central potential 
of Eq.(8) after setting Z — rja'^, B — rfa'^ and C = 0. Therefore, the energy spectrum and 
wavef unctions, for this system can be written immediately from Eq.(32) and Eq.(33) as 

E ~ - -/^(^^'^')' 



2r 

and 

*n,n,m(r,^,9^) = A^nr's-^"^Lf +^n26nr)sinV^^^'^(^)pi/^^^ 

(36) 



respectively, with I — fi -\- ^Jm? -\- (r^cr)^, n — 0, 1, 2, • • • . This result agrees with that of 
in Refs. [34,38] where energy spectrum has been calculated by solving Schrodinger equation 
using the path integral solution. 



2. Aharonov-Bohm plus Dime monopole potential 

We consider the general case involving a vector potential A = '^^e^^v^ ''^^^ -^(^) ~ 
$/27r + g{l — cosO), where $ is the magnetic flux, together with a scalar potential V{r) = 
— This problem of a relativistic Dirac electron in the presence of a Coulombic field plus 
Dirac magnetic monopole and the Aharonov-Bohm (AB) potentials have been discussed by 
Hoang et al. [39]. Both potentials can be chosen to be of arbitrary strength depending on 
the values of the coupling constants $ and g. The effective potential for this system is 
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Vee{r,e) 



+ 



a — b 
sin^ 9 



(37) 



where a = e^F'^{6) and h = 2meF{9). This potential is a special case of the non-central 
potential of Eq.(8) with A — — Ze^, B — a — b, and C — 0. Now if we define 



q — —qe, m — q — m, 

^ 27r 



(38) 



this gives 



B — {171 + q cos 9)'^ — m} 



(39) 



Therefore, Eq.(23) may be rewritten as 



d^Q{e) 2s dQ{e) ^ 1 



ds (1 - s^y 



q^ + /(/ + 1)J - 2qms + /(/ + 1) - m^J e(^) = 0. 

(40) 



By defining the following associated polynomials 

f{s) = -2s, a{s) = 1 - s^ a{s) ^ - + l{l + 1)] - 2qms + l{l + 1) - m^ (41) 

and using Eq.(6), 7r(s) can also be found in four possible values as 

(|m|s+|g|), /ci = Z(Z + 1) + - m^, if < |m| , 

— (|m| s + |g|) , ki = l{l + 1) + q^ — m^, if |g| < |m| , 

{\q\s+\m\), k2 = l{l + l), if \m\ < \q\ , 

-(|g|s + |m|), A;2 = /(/ + !), if |m| < |g| , 

where the roots are determined by means of the same previous procedures. The r(s) is 
obtained as 



7r(s) 



(42) 



t{s) = < 





< 


m 




< 


m| 



2[(|g| - l)s+ |m|], k2^l{l + l), if \m\ < |g| , 

-2[(|g| + l)s+ |m|], A;2 = /(/ + !), if |m| < \q\ . 



(43) 
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Imposing the condition t'{s) < 0, we get two particular solutions: 

i) For the case ki = l{l + 1) + — m^, we obtain A = /(/ + 1) + — — |m| and 
A~ = n{n + 1) + 2n \rn\ , which in turn give the following quantum number: 



1 

"2 + 



n + \m\ + - 



n 2 



q\ n = 0,l,2,--.. 



(44) 



and hence energy spectrum, Eq.(18), with the previous quantum number becomes 



E. 



n,n 



2r 



n+l + Mn+\m\ + iy 



2 ' 



\q\ < \m\ , 



(45) 



and the wavefunction is given in Eq.(28) with the allowed quantum numbers defined by 
Eq.(44). 

ii) For the case k2 = l{l + 1), we obtain A = /(/ + 1) — \q\ and A~ = n{n + 1) + 2n \q\ , 
which in turn give the following quantum number: 



' = -2 + 



n 2 



n+\q\ + 



-q\ n = 0,l,2, 



(46) 



Hence, the energy spectrum, Eq.(18), with the previous quantum number in Eq.(46), be- 
comes 



TP ~ 

n.n 



2r 



i> m > M 



(47) 



Further, the wavefunction is given in Eq.(30) with the allowed quantum numbers defined by 
Eq.(46). 

B. A generalized Oscillatory potential 



The well-known ring-shape oscillator potential 



V{r, 9) = A'r' + 



B 



(48) 



r2 sin^^' 

where B being strictly positive constant. The limiting case C = corresponds to an isotropic 
harmonic oscillator and will serve as a testing for the results to be obtained. 
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The Schrodinger equation given by Eq.(9), for a particle in a ring-shape oscillator po- 
tential can be reduced to an ordinary differential equation in the radial part: 



2 2 K^ + l)^ 



21 



and another one for the angular part: 

d^e{e) cose dQ{e) 



+ - 



sin 9 dO 



+ 



^(^) = 0, g{r)^rR{r), 



(49) 



(m^ + B) 



(50) 



We shall first solve the radial part of the Schrodinger equation. Employing a simple trans- 
formation s — r^, we get 



9"(s) + ^^g'is) + 



-a'^s'^ + £ls-l{l + lj 
4^2 



9{s) = 0, 



(51) 



with 



2 2nEn 2 2/xyl2 



— i-2 



{En > 0), a' = 



Further, defining the following associated polynomials: 

f (s) 1, a{s) = 2s, a{s) = -a^s"^ + els - l{l + 1), 

and using Eq.(6), 7r(s) can also be found in four possible values 

as + 1 + 1, ki^ ^ + a{l + ^), 



(52) 



(53) 



7r(s) = < 



(54) 



-{as + l), ki = ^ + a{l + l), 
as -I, k2^^-a{l + ^), 
-as + 1 + 1, k2^^-a{l + l), 

where ki^2 is determined by means of same procedures as before. The t(s) is also obtained 



as 



t{s) 



3 + 2as + 2l, ki = '-f + a{l+l). 



l-2as-2l, ki^ ^ + a(l + l), 
l + 2as-2l, /e2 = f -a(/ + i), 
3-2as + 2l, k2 = ^ - a{l+ I). 



(55) 



13 



Imposing the condition t'{s) < 0, one gets the following physical solution: For the case 

2 2 

k2 — Y^ '^(^ + i)' '^^ obtain have A„ = 2an and A = ^ — a{l + |). With A — An, the energy 
eigenvalues obtained are [40]: 



nl 



. — (4n + 2/ + 3)^, n = 0,1,2,---. (56) 



Further, to find eigenfunctions, we first determine the weight function which satisfies the 
conditions in Eq.(53) and Eq.(55), using Eq.(5) we get 

p(r) = e-"'-V(''+^\ (0 < r < oo). (57) 

Substituting Eq.(57) into Eq.(4), we obtain the first part of the wavefunction 

yn{r) ^ C^L^-^'/'\ar% (58) 

where C„ being the normahzation factor. Further, using Eq.(7), the other part of the 
wavefunction is 

0„(r) = r'+^e-"'^'/". (59) 
Finally, the radial part of the wavefunction has the general form: 

Rniir) = Nye-^^'/'Lj+'/'\ar'), (60) 

where is a new normalization factor. It is clear from Eq.(60) that this case is a true bound 
state solution [35] . The total wavefunction can be obtained from the combined solutions of 
Eq.(30) after setting C — and also from Eq.(60) as 

where A„ being a new normalization constant. 

IV. RESULTS AND CONCLUSIONS 



In this work, the Schrodinger equation with a class of non-central but sepcrablc potentials 
has been studied. The analytical expressions for the energy levels and eigenfunctions have 
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also been found analytically for every case. This method is simple and useful in solving 
other complicated systems analytically. Therefore, it is found that the energy eigenvalues 
are consistent with the results obtained by using other methods. We also point out that these 
exact results obtained for the non-central potentials may have some interesting applications 
in the study of different quantum mechanical systems and atomic physics. 
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